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Abstract: We consider a model in which dark matter is a composite baryon of a dark 
sector governed by SU(3) gauge theory, with vector-like quarks also charged under £/(l)y. 
The model provides simple answer to the dark matter stability problem: it is a result of the 
accidental dark baryon number conservation. And with an analogy to QCD, all physical 
quantities of the dark matter can be calculated by rescaling the QCD experimental results. 
According to the thermal freeze-out mechanism the mass of the dark matter is predicted 
to be 0(100) TeV in order to achieve a correct relic abundance. Such heavy dark matter 
is in general hard for detection due to small dark matter number density in the universe. 
However, dark baryon number in our model is not necessarily strictly preserved thanks 
to operators suppressed by the Planck scale, and such decay operator results in a decay 
lifetime marginal to the current detection bound. We show our model with 0(1O 27 ) s dark 
matter decay life time can explain the AMS-02 anti-proton data, if it is experimentally 
interpreted as an access, although some theoretical uncertainty may weaken its significance. 
We also investigate other phenomena of this model such as the extragalactic gamma ray 
and neutrino signatures. 
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1 Introduction 

A lot of particle physics candidates have been worked out for dark matter (DM), which 
takes about 26% of the energy density of our universe [1]. All of the models have, at least 
after tuning to certain parameters, to explain some basic facts of DM such as its stability, 
weakness of its interaction with the standard model (SM) matter, its relic abundance and 
so on. Mainly motivated by the stability problem, we revisit the composite DM model, 
in which DM is a baryonic or antibaryonic composite particle in a hidden strong gauge 
interaction, or the technibaryon DM which has a long history [2-4] (see more recent similar 
works [5-10], and more generally the references in [10]). 

In the simplest put our composite DM model is a copy of quantum chromodynamics 
(QCD) at different scale. Some hidden strong gauge group at certain scales has exactly the 
same QCD phenomena of color confinement and chiral symmetry breaking, they contain 
fermionic dark baryons as composite bound states in totally antisymmetric representation. 
In particular we will see that our model also has an SU( 3) strong gauge group, which 
will be introduced below. Protected by the accidental baryon number, the QCD proton is 
stable empirically, even if the theory is promoted to the grand unified theory (GUT) and 
proton decay is allowed. The stability of (the lightest) dark baryon is similarly guaranteed, 
as a consequence of the accidental dark baryon number conservation. 

Within baryonic DM models there are still a lot of possibilities, and some of the 
models are actually difficult to test. If the dark sector and the SM sector are connected 
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only by some portal particle, then the typical interaction between the two sectors may 
be too week to see (e.g. the direct detection cross section). Moreover, for reproducing 
the correct relic abundance we would like to stick to the simplest assumptions that the 
DM is a thermal relic (as well as implicitly the dark baryon antibaryon are symmetric in 
number). However, in the following we will see that the repetition of the QCD exactly 
provides us a tool to calculate the thermal freeze out cross section, and the DM mass is 
therefore determined to be 0(100) TeV [6-10]. Such high scale DM mass implies a very 
small DM number density in the universe, and suppressing any DM event rates such as 
direct detection and annihilation, making it invisible. In particular, we have checked that 
even with the largest possible annihilation cross section of the s wave unitarity bound, the 
most optimistic NFW DM profile and large exposure time assumptions, next generation 
gamma ray indirect detection experiment such as the CTA still cannot probe/constrain 
such a high scale DM candidate from the galaxy center, similar to [11]. 

Still we manage to test our model. Instead of connecting the two sectors with a portal, 
we consider the following criteria [8] that, while the dark baryon itself must carry no net 
SM charges for being sufficiently weakly interacting with the visible sector, the vector-like 
quarks in the hidden sector are charged under some SM gauge groups of SU (3) c x SU (2) l x 
U(1)y- And the second problem of being too heavy on the other hand weakens the DM 
stability if it is not absolute, eventually enhance the decay rate to a detectable level as a 
benefit. For our benchmark 5C/(3)hid x SU(2)r x U(1)b~l model the decay spectrum can 
be determined and the decay rate happens to be marginal to the current bound, which 
eventually makes the model testable. For example we compare our model prediction of 
DM decay against the preliminary AMS-02 antiproton data, as well as checking with the 
Fermi extragalactic gamma ray and the IceCube neutrino data. In all, (unlike many early 
work) we have provided answers to every aspect of DM physics, by building a model which 
is simple, and detectable or falsifiable in the near future. 

This paper is structured as followed. Section 2 is a quick check of the stability in 
more details. In Section 3 we discuss the rescaling calculation of the DM annihilation cross 
section, and eventually show the DM mass scale of 150 TeV. Next Section 4 is devoted to 
the discussion of the specific gauge group choice, in which an 0( 10 27 ) s decay life time is 
estimated and the SU (3)hid x SU( 2)r x U(1)b~l is justified, with discussion of constraints 
(such as the one from the big bang nucleosynthesis) for various scales implicitly in the 
model building. We check the recent AMS-02 antiproton data with the decay life time 
as a major test of the model in Section 5, as well as check other constraints such as the 
extragalactic gamma ray and neutrino flux. At last we conclude in Section 6. 

2 The Dark Matter Stability 

The solution of DM stability, or suppression of any possible decay operators by composite¬ 
ness, still allows a range of hidden gauge groups. Suppose the hidden strong gauge group 
is an SU(N ) gauge group, with N an odd number for baryon to be a fermion. The dark 
baryon is in the totally antisymmetric representation of N dark quarks, and as the proton 
decay in GUT the effective decay operator is at least N + 1 fermions times together, with 


- 2 - 


the new “1” being some (not necessarily to be new, as we will see) “lepton”. To be overall 
mass dimension 4 we need high scale A suppression of power | (N + 1) — 4. In such analysis 
the DM has an decay width of order 

T ~ m 3 x N ~ 4 /A 3N - 5 . (2.1) 

This hidden strong gauge group has no relation with the postulated GUT in the visible 
sector, so the first choice of the large scale is the Planck scale. Then we can see that even 
for a minimal N = 3, the DM can be stable enough. However, the choice of N = 3 will 
actually turn out to be tricky and interesting, as shown in Section 4. 

Further comment is that we expect this dimension-6 QQQL operator for N = 3 gives 
leading contribution to DM decay, simply because given the dimension of 3 x | (with | 
being the quark dimension) which must appear as a whole in any possible decay operator, 
it is unable to construct any operator with even lower dimension. Within the framework of 
dimension-6 QQQL operator the decay process exactly resembles the QCD proton decay, 
and the lattice calculation of the QCD proton decay matrix element can be rescaled to give 
results in the dark sector, given that the lepton L mass are much lower. Specific choice 
such as whether the L is charged or neutral will depend on the model building, which will 
be addressed in Section 4. Apparently going to even higher dimensions there are various 
operator choices such as just adding powers of the SM H to the operator, but they will 
have much more suppressed decay contributions and are irrelevant. 

3 Dark Matter Annihilation and Mass 

Due to the strong dynamics nature there is no way to calculate from the first principle the 
annihilation cross section of baryons with antibaryons at the low velocity regime. A useful 
first approximation is given by the unitarity bound [12] 

4vr(2^ + 1) 

CTUrei ~ -- 

m B Vre 1 

This unitarity bound can be obtained by solving a nonrelativistic Schrodinger equation, 
with a complex potential to decrease the exiting wave amplitude. Achieving the unitarity 
bound corresponds to a vanishing exiting i partial wave, which is a result of alignment 
of the potential configuration. Potential scattering is nonperturbative in nature, in the 
Feynman diagram point of view messenger particles are repeatedly exchanged, forming 
ladder diagrams. And enhancement of the hard scattering/annihilation process are taken 
into account in this way, recovering the Sommerfeld enhancement [13, 14]. 

The approximation of Eq. (3.1) is indeed supported by the real QCD experiments. In 
the low velocity regime higher partial waves are suppressed by powers of velocity, leading 
to a general s wave domination if it is not forbidden by symmetries. However for a thermal 
freeze out velocity of n re i — 0.4, the s wave may be further assisted by a p wave or even 
higher partial waves. In Fig. 1 we see from experimental results that the annihilation cross 
section is indeed higher than s wave unitarity bound by a factor of a few, due to mainly 
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Figure 1. The QCD measurement of baryon antibaryon annihilation cross section av in the low 
relative velocity region. At higher relative velocity av will exceed the s wave or the (s + p ) wave 
unitarity bound, pp data are taken from [15-18], and rip data are taken from [19, 20]. 

the p wave contributions. A significant p wave contribution same to our case has already 
been shown in [20]. 

In Eq. (3.1) the true QCD annihilation cross section scales as j 2 with tub the QCD 
baryon mass, while the 2i + 1 adopt an effective value of about two to four (note that the 
annihilation remains rather constant with variation of relativity velocity in Eq. (3.1), which 
means the effective 2^ + 1 has a rough scaling with u re i, to cancel the u re i dependence). If 
it is applied to the dark sector baryon antibaryon annihilation the same effective 2^+1 
should be kept, difference is only in replacing the QCD baryon mass (about 1 GeV) by the 
dark baryon/antibaryon mass. Namely the DM annihilation cross section is just a rescaling 
of the QCD annihilation cross section 

° V re\( m x) = ™rel X (““) ’ ( 3 ‘ 2 ) 

where is the averaged value of purple and blue points in the Fig. 1. In the following 

we will see that our benchmark model is also based on the SU(3) strong gauge group; and 
we are taking implicitly that the current quark mass is at least one order smaller than 
the baryon mass, making the chiral symmetry breaking effect small as the QCD. In such 
case the effective theory can be approximated by only one scale, which is the QCD scale 
of typical baryon mass for the real QCD, and the dark baryon mass for the dark sector 1 . 
Thus the strong dynamics of the dark sector is quite similar to the real QCD case, and we 
expect the rescaling calculation to be very good. 

1 The other scales introduced in Section 4 . 2 , 4.3 has nothing to do with the annihilation. 
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The thermal freeze out mechanism is the most elegant way to reproduce the observed 
relic abundance, in which we only set the DM annihilation cross section to the famous 

3 x 10 -26 cm 3 s _1 , corresponding to fl x h 2 ~ 0.12. When we do this to the baryonic dark 
sector we get 

m x ~ 150 TeV. (3.3) 

Here we make a few comments and comparisons on the result, which is crucial to the 
following analysis. This number is larger than the literature value of the s wave unitarity 
bound of the DM mass in thermal freeze out, which calculated from Eq. (12) of [12] is 
about 86 TeV for a Dirac fermion. Note that the latter number is assuming a freeze out 
relative velocity of about 0.4c. In Fig. 1 at such relative velocity we can read the true 
baryonic annihilation cross section is almost four times larger than the s wave unitarity 
bound, so a further factor of two difference of DM mass can be understood. The discussion 
for higher partial waves in [12] is implicitly assuming that one partial wave will dominate 
the total annihilation cross section. This assumption fails for a strong dynamics system 
such as the QCD case, which has already been shown in Fig. 1 for the real QCD. Still in 
the same small relative velocity regime due to the strong dynamics nature many partial 
waves contribute to the annihilation amplitude and possibly interfere with each other, in 
that case the implicitly omitted 2^ + 1 factor in [12] makes a difference. 

4 The Benchmark SU (3)hid x SU( 2)# x U(1)b-l Dark Matter Model 

4.1 The Hidden Gauge Group 

The high scale of DM does provide one advantage, that the increase of decay width from 
Eq. (2.1) makes better chance of its detection in the decay channel. We start with checking 
N = 3, a reduced Planck scale 2 square suppression and the dark baryon mass of 150 TeV. 
Eq. (2.1) gives T ~ m 3 / Apj ~ (150 x 10 3 ) 5 /(2.4 x 10 18 ) 4 = 2.3 x 10 -48 GeV, which 
corresponds to a decay life time of 2.9 x 10 23 s. This estimation has already revealed an 
amazing coincidence, that the minimal model leads to a decay life time close to the current 
observational limit. 

We can use lattice calculation of the GUT proton decay to further improve the esti¬ 
mation. Rescaling Eq. (10) of [21] of such calculation gives further significant corrections 
of an extra factor of 1/(327 t) from the phase space, and another 10~ 2 or so as the lattice 
calculated matrix element amplitude square of — (0.1 GeV 2 ) 2 divided by the QCD 
proton mass to the fourth. Putting altogether the decay life time is a few times 10 2 ' s. 

The choice of hidden gauge group of 517(3) has additional benefit, of further justifying 
the previous rescaling of the annihilation calculation. The real QCD has strong gauge 
group of SU(3) C , or baryon consisting three quarks. The choice of 5f7(3)hid also preserves 
this structure, which is expected to minimize the uncertainty of our ignorance of strong 

2 Here we use the reduced Planck scale ( 8 - 7 rG ) _ 2 = 2.4 x 10 18 GeV rather than the Planck scale G “ 5 = 

1.2 x 10 19 GeV. The reduced Planck scale appears naturally in supergravity theory of particle physics. The 
Planck scale, on the other hand, already gives too low a decay rate. 
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dynamics, making the rescaling the most reliable compared with other choices of the strong 
gauge group. 

4.2 Introduction of the SU(2)r x U(1)b-l 

Let’s take a look at other aspects of the minimal SU( 3)hid model. The requirement that it 
is also charged under some SM gauge group is minimally satisfied for the SM U(l)y group, 
and that will induce the SM electromagnetic interactions and the current quarks will be SM 
electromagnetically charged. To further build our model we refer to the most common way 
of extension of the U( l)y in beyond SM model building, namely the SU(2)r x U(1)b-l ~► 
U(1)y left-right symmetry model, which also gives the neutrino masses. An SU(2)r gauge 
group behaves as the counterpart of the SU(2)l isospin, and the U( 1) before this breaking 
becomes B — L. The hidden quarks are postulated to be in the doublet representation of 
SU( 2) r . 



Gauge 

SU(3) h 

SU( 3) c 

SU( 2) R 

SU(2) l 

U(1) B -l 


$ 

1 

1 

3 

1 

+1 

New Particle 

Ql 

3 

1 

2 

1 

+5 


Qr 

3 

1 

2 

1 

+* 


H 

1 

1 

2 

2 

0 


ql 

1 

3 

1 

2 

+ 5 

SM 

QR 

1 

3 

2 

1 

+ 5 


II 

1 

1 

1 

2 

1 

2 


Ir 

1 

1 

2 

1 

1 

2 


Table 1. The particle content and their quantum numbers. 


Going back to the minimal dimension-6 decay operator QrQrQrIr, the “lepton” is 
also in the doublet representation of SU(2)r. Actually it can be identified as the SM right 
handed charged lepton, which is also promoted by the beyond SM SU(2)r. The model 
is shown in Table 1. Here we also introduce an SU(2)r triplet 4>, which gets some large 
vacuum expectation value to break the SU(2)r , and gives large masses to the SU(2)r 
W boson and the right handed neutrino [22]. The right handed neutrino mass scale is 
postulated to be determined by the natural seesaw mechanism, or much larger than the 
DM scale of 150 TeV, so the Ir could only be the charged right handed lepton, in order not 
to be kinetically forbidden. We do not expect forbidden of right handed neutrino decay 
channels significantly affects the estimation of the DM life time. 

At the level of the elementary particle our model can be summarized as the renormal- 
izable Lagrangian 


C D £lr — (A t^e^lR + h.c .) 

_ I<7 a n a i iv — \w a W a ^ u --F F^ v 
+ Q l (0 - M)Q l + Q R (iIp - M)Qr + {D^D^. 


(4.1) 
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The £lr is the left-right symmetric extension of the SM, and we write the allowed neutrino 
Majorana mass terms separately. Note the e = ia 2 = (_i 4 ) in the SJJ(2) contraction, and 
in this matrix basis (<f>) = (y q). Each covariant derivative can be read from Table 1. 

Note that if we choose the SM SU(2)l rather than the SU(2)r for the dark quark, 
then the dark baryon is also charged under the SM SU(2)l. This still suits the definition 
of “WIMP”, but the scattering with ordinary matter in direct detection will have charge 
induced W 2 * 1 scattering rather than the current dipole induced photon scattering as the 
leading one, and we do not consider this situation in our scenario. 


4.3 The Mass Splitting between Dark Neutron and Proton 

After the introduction of SU(2)r, the baryonic DM candidate is also promoted into doublet. 
The dark neutron is the DM candidate. The dark proton need to be heavier than the dark 
neutron, and after thermal freeze out decay early enough to it. 

Dimensional analysis suggests such decay has a width of T oc Am 5 /my/ [23], where 
mw is the SM (or SU(2)r) W boson mass for real QCD (or the dark sector), and Am is 
the mass difference between the SM (or dark) proton and neutron. The simple requirement 
that the dark proton to dark neutron decay life time is before the big bang nucleosynthesis 
(< 1 s) corresponds to a constraint of 


Am > 0.22 


( m w R \ t 
V 10 s GeV 7 


TeY, 


(4.2) 


which is obtained by rescaling the QCD free neutron to proton decay. 

The mass difference between the dark proton and the dark neutron comes in two ways: 
One is the electromagnetic radiative self energy correction which applies only to charged 
proton. The other is the current quark mass differences. In real QCD the two happen to 
cancel with each other, with the latter the dominant one which makes neutron heavier. 
However in our dark baryonic sector the two can be additive. The radiative self energy 
correction will contribute a fixed mass splitting of about 0.1% of the baryon mass, which 
alone will be sufficient for a light m\y R < 10 8 GeV or a low seesaw scale, nothing else is 
needed. 

On the other hand, for a heavy mw R the radiative self energy correction is not enough, 
and the current quark mass splitting shall be the dominant source of mass differences 
between dark neutron and proton. In our setup it can be generated by the dimension- 
5 operator Q&&Q, in which the vacuum expectation value of indeed makes the “up 
type” dark quark as well as the dark neutron heavier (viz. we give a theory for the up- 
down current quark mass splitting which is unexplained in QCD). Assuming the SU(2)r 
coupling is order one, the vacuum expectation value of is expected to be at the same 
order of m\y R - Here we need to introduce a new effective suppression scale A' for the 
dimension-5 operator. In Fig. 2 we plot the mw R vs. Am with several A's. We can see 
that the big bang nucleosynthesis constraint only exclude some region of A' > 10 14 GeV, 
even when the Am is reasonably small. 
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Figure 2. The big bang nucleosynthesis constraint of r < 1 s (shaded region excluded), on the dark 
neutron proton mass splitting Ato = 0.15 TeV + /A' (electromagnetic radiative self energy 

correction plus current quark mass splitting from dimension-5 operator) with mw R ■ 

4.4 The Dark QCD Spectrum 

In addition to the baryonic new particles in the dark sector possibly there are also mesonic 
new particles. Among them the most interesting ones will be the lightest, which are 
the pseudo Nambu-Goldstone (pNG) bosons 3 . Corresponding to the QCD pions in the 
hypothetical proton decay, they can be the other decay product of the dark neutron except 
for the right hand charged leptons. Like the SM pions they are not stable either, the 
charged one can decay through a virtual Wr into SM fermions, with the chirality flipping 
mechanism making the top bottom pair the dominant channel, and the neutral one will 
dominantly decay into photon pairs through anomalous triangle diagram. 

In real QCD the pion mass is determined by the Gell-Mann-Oakes-Renner relation [24] 
from the partially conserved axial current principle 

2 (m u + m d ) , 

= -72-w<?) 

J IT 

In our dark QCD with the baryonic particle scale of 150 TeV, we expect the — {qq)/f% 
also at the scale of hundred TeVs. The current quark mass is a complete free parameter, 
making the dark pNG pion mass also a free parameter. On the experimental side, the 
charged component of the dark pion should decay through a virtual Wr which is exactly 
like the QCD pion decay through a virtual SM W, and chirality flipping mechanism makes 
the tRbR channel the dominant one. However the usual searches for W'± —>• tb (or the 
H^ —>• tb which shares the same final state) assumes a singly produced n^, and in our 
model the dark pion should be produced in pairs by a very energetic gamma. Since searches 
for the top and bottom dijets in pairs have not been performed in this topology, and the 

3 We assume some other broken symmetry to protect this Nambu-Goldstone bosons from having large 
masses. Note that one cannot equate the SU(2)r as the isospin SU(2)l in real QCD the chiral breaking 
of which leads to the QCD pion, for the quarks in our scaled-up QCD are vector-like. 


(4.3) 
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production is through electroweak process which is very suppressed compared to a QCD 
process, there is essentially no experimental constraint. 

We also expect other mesons in the spectrum which differ in discrete quantum num¬ 
bers, such as the dark p meson, the p meson, the u meson and so on. Without a universal 
understanding of the spectrum in strong dynamics we do not know the exact mass com¬ 
parison with the dark neutron, so we will ignore their contribution in the dark neutron 
decay. 

5 DM Decay Detectibility 

Next we will compare the theoretical prediction of the 0(1O 27 s) DM decay life time to 
current and future experiments, to see whether they are detectable. We find that the AMS- 
02 preliminary antiproton data indeed favors an O(10 27 s ) DM decay life time, and it does 
not exceed the Fermi EGB bound and the IceCube neutrino bound. As a benchmark we 
ignore all the SM particle mass including the top quark which are much smaller, and choose 
the dark pion mass to be 1/10 of the dark baryon mass or 15 TeV. We want to emphasize 
that it should not be taken as a rigorous data driven fit. 

Except for the three DM indirect detections, we found direct detection experiments 
are also constraining and complementary. Our result agrees with [10] and we would like to 
redirect reader to their paper but not repeat it here. 

5.1 Proton and Antiproton ratio 

The dark neutron will decay through the proton decay like chain of N —>• 11“ — > 

as discussed before. The neutral decay mode of replacing the by the (on-shell) right 
hand antineutrino is kinematically forbidden. The branching ratio for i + is assumed equal 
for three families. The same but every particle replaced by antiparticle decay chain applies 
to dark antineutron. 

For the most interesting decay product of top quark, the energy distribution can be 
determined analytically in the sequential two body decay. With mn = 15 TeV fixed the 
charged lepton has a fixed energy of Eg = (m 2 — m 2 1 )/(2m x ) = 74.25 TeV, and the top and 

bottom will be evenly distributed in the energy region from E qm j n = Eu — E ^ — m u) 

to E qmax = \{Eu + E ^ — mpj), where Eu = (m 2 + m^ l )/(2m x ) = 75.75 TeV. All the 
prompt decay spectra are then calculated by 



where the primed are taken from the PPPC4 [25] for i = pp, e + e ,7 and so on. 4 

4 The PPPC4 cosmic ray spectrum calculation is only available for unpolarized primary quarks, but in 
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Figure 3. The antiproton to proton ratio fitting to the AMS-02 preliminary data for the MAX (left 
panel) and MED (right panel) propagation parameters. The DM alone, background alone, and sum 
of DM and background contributions are shown as green dashed-dotted line, blue dashed line and 
the black line respectively. 


The antiproton prediction are shown in Fig. 3, to the AMS-02 preliminary antiproton 
data [26] which shows a slight excess but still can be viewed as consistent with back¬ 
ground. In this work we use the the NFW halo profile and the MAX/MED propagation 
parameters [27, 28] for antiproton, as argued recently in [29] that the MED is less favored by 
isotopic abundance data from AMS-02, but for completeness we still leave our prediction 
for MED on the right panel. The background is taken from [30]. Considering K > 10 GeV, 
the best-fit DM decay time is r = 2.4 x 10 2 ' second for MAX and r = 4.5 x 10 26 second 
for MED. The chi-square improvement of (x'g um — Asecondary)/^- 0 -/- can achieve 19.97/17 
(66.29/17) for MAX (MED). 

We have checked that varying the dark pion mass (e.g., by a factor of two) has negligible 
effect on the spectrum (Eq. 5.1), as well as the experimental statistical strength. The result 
always holds as long as the dark pion mass is much smaller than the dark matter mass. 

5.2 Extragalactic Gamma Ray and Neutrino 

The EGB :) prediction comparing with the Fermi-LAT 4 year data [33] are shown in Fig. 4. 
We ignored the local DM decay contribution such as that from our Milky Way, only working 
for redshift z m j n = 10 -4 to z max = 2 where the gamma ray is effectively cut off by the optical 

our model the primary tn and bn are all right-handed. In parton shower and hadronization process the 
helicity induced difference in yield spectrum is arguably small (which can affect a few percent only in the 
high energy region), hence we directly use the unpolarized result. 

5 In [31] it is pointed out that the most stringent gamma ray limit for decaying DM is from the angular 
cross-correlation of low-redshift sources, but the improvement, c.f. [32], seems weaker in our high energy 
region. 
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Figure 4. The EGB gamma ray (left panel) and neutrino (right panel) energy spectrum predictions 
based on two benchmark decay times. The gray and green bands represent the 2a favored decay 
time by AMS-02 antiproton data with propagation MED and MAX respectively. 


depth for all interested energy. We also ignore the contribution from inverse Compton 
scattering of our charged decay product with the CMB photons, which is less important in 
our interested high energy region and makes the result conservative. The gamma ray flux 
can be described by 


= r«n~ e -T(z,E) d N 

dE 4vr m x r x ,A m . n “ H(z) dE 1 ’ ’ 


(5.2) 


where H(z) — (H x + Hft)(l + z) 3 + Ha is the Hubble function and the data of Hq, fi x , 
and critical density p c are taken from [34]. The redshifted photon energy E measured at 
earth is related to the initial energy of E z = E(1 + z) at production, and the optical depth 
suppression e~ r ^ z,E ^ is taken of the min UV case in PPPC4 [25]. For the background we use a 
power law with cutoff, ignoring possible astrophysical sources such as blazars, star-forming 
galaxies [35], and misaligned active galactic nuclei [36, 37]. 

In general we find an 0( 10 27 ) s decay life time is consistent with the Fermi EGB 
gamma ray data. For example, with r = 2.4 x 10 2 ' s favored by the AMS-02 antiproton 
data with MAX propagation, our model prediction is well below the Fermi EGB gamma 
ray data. On the other hand, the r = 4.5 x 10 26 s favored by antiproton data with MED 
propagation is almost excluded, but still barely allowed if including the 2 cr error band of 
the antiproton data. 

Also the 0(1O 27 ) s decay life time is allowed in the neutrino channel. There is no recent 
isotropic neutrino data in our desired energy regime 6 , so we use an extrapolation [38] fitted 

6 In our model we cannot have primary SM neutrinos and all produced SM neutrinos are secondary, and 
their energy should be much smaller than the DM scale of 150 TeV, so unable to account the IceCube 100 
TeV to PeV neutrino signals. 


, Ha 
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from the IceCube neutrino data at an even higher energy spectrum [39], assuming a single 
unbroken power law and equal neutrino fluxes of all flavors. In the right panel of Fig. 4, it 
is shown that the 0( 10 27 ) s decay life time is at least two orders below the current IceCube 
neutrino extrapolated bound. 

6 Summary 

We have proposed a new scenario of baryonic DM based on a strong hidden SU (3) gauge 
group, which also connects to the visible sector through the SU (2)r x 17(1) b~l gauge group. 
This composite gauge group is primarily motivated by the DM stability consideration, but 
the hypothetical Planck scale suppressed decay can induce interesting possible signals while 
the connection sector SU(2)r x U(1)b~l is added. 

In the minimal model building we introduce the following scales 

• The DM or the dark baryon mass scale of 150 TeV, which is determined by rescaling 
the real QCD baryon antibaryon annihilation cross section to the “WIMP miracle” 
cross section. 

• The dark pion mass scale which is much lower. 

• The SU(2)r breaking scale, which is free but appears in the natural seesaw mecha¬ 
nism. 

• The optional free effective suppression scale of the dimension-5 operator. 

We find that theoretical estimation and fitting to recent preliminary AMS-02 antipro¬ 
ton data are both consistent with a decay life time of a few 10 2 ' s. Although being heavy, 
our DM model has a perspective of detection in the near future. 
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